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The  terminal  guidance  problem  for  missiles  intercepting  maneuvering  targets  with  terminal  impact  angle 
constraints  is  investigated.  Regarding  the  target  acceleration  as  an  unknown  bounded  disturbance,  novel 
guidance  laws  based  on  integral  sliding  mode  control  (ISMC)  method  technique  are  developed.  The 
first  one  is  a  linear  integral  sliding  mode  (ISM)  guidance  law,  which  can  guarantee  the  line-of-sight 
(LOS)  angular  rate  and  the  LOS  angle  asymptotical  convergence  with  infinite  time.  To  further  improve 
the  convergence  characteristics  of  guidance  system,  a  nonlinear  ISM  guidance  law  is  developed,  which 
guarantees  the  LOS  angular  rate  and  LOS  angle  finite-time  convergence  characteristics.  However,  to 
guarantee  the  guidance  system  has  a  good  performance  for  dealing  with  target  acceleration,  the  switch 
gains  of  both  linear  and  nonlinear  ISM  guidance  laws  need  to  be  chosen  larger  than  the  bound  of  the 
target  acceleration.  It  will  lead  to  chattering  problem.  To  reduce  the  chattering  phenomenon  and  improve 
the  performance  of  system,  nonlinear  disturbance  observer  (NDOB)  is  employed  to  estimate  the  target 
acceleration.  The  estimated  acceleration  is  used  to  compensate  to  actual  target  acceleration.  Then,  two 
novel  composite  guidance  laws  combining  linear  and  nonlinear  ISM  guidance  laws  with  feedforward 
compensation  terms  based  on  NDOB  are  developed.  Finally,  simulation  comparison  results  are  provided 
to  demonstrate  the  effectiveness  of  the  proposed  methods. 
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1.  Introduction 

The  objective  of  a  missile  guidance  law  is  to  provide  accelera¬ 
tion  commands  to  the  missile  autopilot  that  produce  a  minimum 
miss  distance  with  respect  to  a  given  target  [32],  For  the  past 
decades,  proportional  navigation  (PN)  guidance  law  is  the  most 
popular  and  widely  used  guidance  law  in  the  practical  missile  in¬ 
terception  engagements  because  of  its  ease  of  implementation  and 
high  efficiency  [19],  The  PN  guidance  law  is  first  introduced  by  [44] 
where  the  target  moves  on  a  straight  course  and  a  plane  curved 
course  are  discussed,  respectively.  The  principle  of  PN  guidance 
law  is  that  the  commanded  missile  normal  acceleration  is  propor¬ 
tional  to  the  line-of-sight  (LOS)  angular  rate  [4],  It  is  well  known 
that  PN  guidance  law  is  the  optimal  guidance  law  with  a  naviga¬ 
tion  ratio  N  =  3,  and  it  can  result  in  successful  interceptions  under 
a  wide  range  of  engagement  conditions  [42],  Then  many  kinds  of 
modified  PN  guidance  laws  are  developed  in  the  literature  [12,24, 
43],  However,  the  PN  and  modified  PN  guidance  laws  are  more  ap¬ 
plicable  for  the  task  of  intercepting  a  nonmaneuvering  target  or 
a  weakly  maneuvering  target.  In  practice,  target  acceleration  can 
change  rapidly.  For  intercepting  a  target  with  powerful  maneuver¬ 
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ing  capability,  a  significant  miss  distance  may  be  resulted  by  PN  or 
modified  PN  guidance  laws  [46],  To  obtain  a  small  miss  distance, 
the  guidance  laws  based  on  modern  control  theory  are  proposed, 
such  as,  optimal  guidance  law  [37],  differential  game  guidance  law 
[26],  nonlinear  Hoo  guidance  law  [38],  Lyapunov-based  guidance 
law  [17],  L2  gain-based  guidance  law  [48],  adaptive  guidance  law 
[10],  sliding  mode  guidance  law,  etc. 

It  is  well  known  that  sliding  mode  control  (SMC)  is  an  effi¬ 
cient  control  method,  not  only  because  it  has  ability  to  deal  with 
nonlinear  perturbations,  external  disturbances  and  parametric  un¬ 
certainties,  which  make  the  system  have  a  good  robust  and  distur¬ 
bance  rejection  performance,  but  also  because  the  discontinuous 
terms  may  allow  the  system  states  to  have  more  quickly  conver¬ 
gence  characteristic  than  continue  control  laws  would  do  [34],  In 
the  guidance  law  design  problem,  the  LOS  angular  rate  is  com¬ 
monly  used  to  define  the  sliding  mode  manifold  because  the  phys¬ 
ical  principle  is  to  attain  a  collision  triangle  [27],  The  zero-effort 
miss  distance  is  another  candidate  for  designing  the  sliding  mode 
manifold  [45],  Considering  the  target  acceleration  as  a  bounded 
uncertainty,  Moon  et  al.  proposed  a  sliding  mode  guidance  law 
[22],  The  target  acceleration  is  suppressed  by  the  target  acceler¬ 
ation  bound.  Shima  et  al.  derived  a  sliding  mode  guidance  law  for 
an  integrated  missile  autopilot  and  guidance  law  design  problem 
with  the  zero-effort  miss  distance  as  sliding  mode  manifold  [27], 
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Lin  et  al.  proposed  an  adaptive  fuzzy  sliding  mode  guidance  law  to 
force  the  missile  to  move  along  the  instantaneous  line  of  sight  [20], 
By  the  fuzzy  rules  learning  online,  the  guidance  law  adjusts  the 
parameters  of  the  fuzzy  rules  and  the  target  acceleration  bound. 
Considering  the  flight  dynamics  of  missile  and  actuator  time  delay, 
Brierley  et  al.  proposed  a  sliding  mode  guidance  law  with  the  dif¬ 
ference  between  the  LOS  angular  rate  and  steering  fin  angle  as  the 
sliding  mode  manifold  [5],  Zhou  et  al.  proposed  an  adaptive  slid¬ 
ing  mode  guidance  law  with  the  product  of  relative  distance  and 
LOS  angular  rate  as  the  sliding  mode  manifold  [47],  Ebrahimi  et  al. 
proposed  an  optimal  sliding  mode  guidance  law  with  the  missile 
terminal  velocity  constraint  by  selecting  the  control  energy  as  the 
optimization  criterion  and  the  combination  of  zero-effort  miss  dis¬ 
tance  and  zero-effort  velocity  as  the  sliding  mode  manifold  [11], 
Alexander  et  al.  proposed  a  sliding  mode  control  guidance  law  for 
integrated  guidance  and  control  system  using  a  sliding  mode  ob¬ 
servers  to  estimate  the  target  states  [1],  All  above  SMC  guidance 
laws  have  the  ability  to  deal  with  the  target  acceleration. 

The  conventional  SMC  is  powerful  in  controlling  nonlinear  sys¬ 
tems,  but,  when  the  system  states  have  not  yet  reached  the  sliding 
mode  manifold,  it  is  sensitive  to  parameter  variation  and  external 
disturbance  [7,35].  In  order  to  improve  the  robustness  in  reach¬ 
ing  phase,  integral  sliding  mode  control  (ISMC)  was  proposed  by 
Utkin  et  al.  [35],  Compared  with  the  conventional  SMC,  the  integral 
sliding  mode  (ISM)  dynamic  starts  from  the  initial  time  instance. 
Therefore,  the  reaching  phase  is  eliminated  and  robustness  of  the 
system  can  be  guaranteed  throughout  the  entire  response  of  the 
system.  Hence,  the  performance  of  conventional  SMC  can  be  im¬ 
proved  by  introducing  ISMC.  Till  now,  ISMC  schemes  for  linear  and 
nonlinear  systems  have  been  studied,  e.g.,  the  design  methods  for 
ISM  manifold  [6,7],  higher-order  ISMC  methods  [16,49],  and  ap¬ 
plied  in  many  control  fields,  e.g.,  induction  motor  [2],  five-degree- 
of-freedom  active  magnetic  bearing  system  [21],  omnidirectional 
mobile  robot  [15],  etc.  The  integral  action  of  ISMC  can  improve 
the  reaching  phase,  at  the  same  time,  it  may  bring  overshoot  and 
lengthen  the  settling  time.  To  overcome  this  problem,  finite-time 
control  strategy  and  NDOB  technique  are  employed  in  this  paper. 

Actually,  the  missile  guidance  problem  considered  in  most  liter¬ 
atures  are  solved  by  the  asymptotic  stability  analysis  which  implies 
that  the  LOS  angular  rate  converges  to  zero  with  infinite  settling 
time.  It  is  well  known  that  finite-time  stabilization  of  dynamical 
systems  may  give  rise  to  a  better  disturbance  attenuation  besides 
fast  convergence  to  the  origin  [3],  Finite-time  control  strategy  for 
the  guidance  systems  have  also  been  studied.  SMC  guidance  laws 
with  finite-time  convergence  were  proposed  by  Zhou  et  al.  [46], 
The  proposed  guidance  scheme  involved  a  continuous  homoge¬ 
neous  function  to  ensure  finite-time  convergence  of  the  LOS  an¬ 
gular  rate.  By  considering  the  dynamics  of  a  missile’s  autopilot  as 
a  first-order  lag,  Sun  et  al.  proposed  a  SMC  guidance  law  with  fi¬ 
nite  time  convergence  [33].  Shtessel  et  al.  proposed  a  composite 
SMC  guidance  law  with  finite  time  convergence  characteristics  for 
a  missile  guidance  system  against  target  performing  evasive  ma¬ 
neuvers  [29]  and  for  a  two-loop  integration  of  guidance  and  flight 
control  systems  with  dynamic  uncertainty  [30], 

For  the  SMC  problem,  there  exists  an  unavoidable  application 
problem-chattering.  To  alleviate  chattering  phenomenon,  one  of 
the  solution  is  to  employ  a  saturation  function  instead  of  signum 
function  in  the  control  input  [36],  However,  to  do  this,  the  distur¬ 
bance  rejection  performance  is  sacrificed  to  some  extent.  Another 
efficient  method  for  alleviating  the  chattering  problem  is  to  em¬ 
ploy  a  disturbance  observer  (DOB)  to  estimate  the  disturbances. 
The  disturbance  estimation  is  used  for  compensation.  The  distur¬ 
bance  observer  technique  is  first  proposed  in  [23],  and  up  to  now, 
DOB-based  control  schemes  for  linear  and  nonlinear  systems  have 
been  studied  and  applied  in  many  control  fields,  such  as  non¬ 
smooth  observer  [18],  robotic  systems  [9],  grinding  systems  [39], 


nonlinear  magnetic  leviation  suspension  system  [40],  mismatched 
uncertainties  [41],  general  systems  [8],  etc. 

When  using  a  missile  to  intercept  a  target,  the  impact  geome¬ 
try  is  often  important.  The  lethality  of  warheads  can  be  increased 
at  a  specific  impact  angle.  Thus,  the  perfect  guidance  law  not  only 
can  guarantee  the  missile  shortens  the  miss  distance,  but  also  can 
guarantee  the  missile  has  an  appropriate  attitude  to  hit  target.  Shin 
et  al.  proposed  a  guidance  law  by  solving  the  two-point  bound¬ 
ary  value  problem  with  impact  angle  constraint  [28],  Ryoo  et  al. 
proposed  energy  optimal  guidance  laws  with  impact  angle  con¬ 
straint  for  arbitrary  missile  dynamics  [25],  Song  et  al.  proposed  a 
time  optimal  guidance  law  in  the  case  of  constrained  missile  ma¬ 
neuverability  and  a  first-order-lag  missile  autopilot  [31],  Harl  et  al. 
proposed  a  second-order  sliding  mode  guidance  law  with  terminal 
constraints  of  impact  angle  and  impact  time  using  a  backstepping 
concept  with  the  difference  between  LOS  angle  and  desired  LOS 
angle  as  sliding  mode  manifold  [13],  Ratnoo  et  al.  proposed  a  PN 
guidance  law  against  a  nonstationary  nonmaneuvering  target  with 
impact  angle  constraint  [24], 

In  this  paper,  a  more  interesting  second-order  nonlinear  missile 
guidance  problem  with  impact  angle  constraint  will  be  considered, 
and  novel  guidance  laws  based  on  ISMC  method  and  NDOB  tech¬ 
nique  are  proposed.  A  linear  ISMC  method  is  employed  for  asymp¬ 
totically  stable  strategy,  and,  a  nonlinear  ISMC  method  is  employed 
for  finite-time  stable  strategy.  The  NDOB  is  used  for  estimating  the 
target  acceleration.  The  proposed  guidance  laws  is  designed  with 
no  information  on  the  target  acceleration,  and,  they  can  yield  ex¬ 
cellent  performance  and  desirable  robustness  properties.  First  of 
all,  a  linear  ISM  guidance  law  based  on  a  linear  ISM  manifold 
inspired  by  Zong  [49]  is  designed  by  using  Lyapunov  theory  un¬ 
der  the  assumption  that  the  target  acceleration  is  bounded  by  a 
constant  upper  bound.  The  law  guarantees  the  guidance  system 
asymptotically  stable.  To  further  improve  the  convergence  charac¬ 
teristics  of  guidance  system,  a  nonlinear  ISM  guidance  law  based 
on  a  nonlinear  ISM  manifold  is  proposed  by  using  finite-time 
Lyapunov  theory.  The  nonlinear  law  guarantees  the  guidance  sys¬ 
tem  finite-time  stable.  However,  to  guarantee  the  guidance  system 
has  a  good  performance  for  dealing  with  target  acceleration,  the 
switching  gain  needs  to  be  chosen  larger  than  the  bound  of  the 
target  acceleration.  It  will  lead  to  chattering  problem.  Then,  to  re¬ 
duce  the  chattering  phenomenon  and  improve  the  performance 
of  system,  an  efficient  disturbance  estimation  technique,  NDOB, 
is  employed  to  estimate  the  target  acceleration.  After  disturbance 
compensation  based  on  NDOB,  the  switching  gain  only  needs  to 
be  larger  than  the  bound  of  the  disturbance  compensation  error 
which  is  usually  much  smaller  than  that  of  the  target  acceleration. 
Hence,  the  chattering  will  be  reduced  while  the  disturbance  rejec¬ 
tion  performance  of  the  guidance  system  can  be  maintained.  Then, 
novel  composite  guidance  laws  combination  of  ISMC  method  and 
NDOB  technique  are  developed. 

The  rest  of  this  paper  is  organized  as  follows.  In  Section  2, 
some  preliminaries  are  briefly  outlined.  The  main  results  are  pre¬ 
sented  in  Sections  3  and  4,  respectively.  In  Section  3,  linear  and 
nonlinear  ISM  guidance  laws  are  proposed.  In  Section  4,  two  com¬ 
posite  guidance  laws  are  developed  combining  the  ISMC  method 
with  NDOB  technique.  To  verify  the  effectiveness  of  the  proposed 
guidance  laws,  simulation  results  are  provided  in  Section  5.  Finally, 
conclusions  are  drawn  in  Section  6. 

2.  Problem  formulation 

This  section  presents  the  mathematic  model  derivation  of  the 
guidance  system  for  missile  intercepting  target.  The  geometry  of 
planar  interception  is  depicted  in  Fig.  1.  We  denote  the  missile 
and  target  by  subscripts  m  and  t,  respectively.  For  the  purpose  of 
guidance  law  design,  missile  and  target  are  assumed  to  be  point 


Z.  Zhang  et  al.  /  Aerospace  Science  and  Technology  31  (2013)  30-41 


Fig.  1.  Planar  interception  geometry. 


masses,  and  only  kinematics  are  considered.  The  corresponding 
equations  of  guidance  system  [26]  are  given  by 


where  f(x)  :  Lf0  ->■  Rn  is  continuous  on  an  open  neighborhood  Uq 
of  the  origin  x  =  0.  The  equilibrium  x  =  0  of  the  system  is  (locally) 
finite-time  stable  if  the  following  two  conditions  are  satisfied:  (i)  It 
is  asymptotically  stable,  in  U,  an  open  neighborhood  of  the  ori¬ 
gin  x  —  0  with  U  c  Uo  C  Rn.  (ii)  It  is  finite-time  convergent  in 
U,  that  is,  for  any  initial  condition  Xo  £  U\{0),  there  is  a  settling 
time  T  >  0,  dependent  on  xo,  such  that  every  solution  x(t,x o)  of 
system  (7)  is  defined  with  x(t,x o)  £  U\{0)  for  t  e  [0,  T)  and  satis¬ 
fies  limt_>T(Xo)  x(t,  Xo)  =  0  and  x(t,  Xo)  =  0,  if  t^T (xo).  Moreover, 
if  U  =  Rn,  the  origin  x  =  0  is  globally  finite-time  stable. 

The  following  lemma  will  be  used  in  the  subsequent  guidance 
laws  development  and  analysis. 


rq  =  \7msin(q  -  9m)  -  Vt  sin(q  -  6t),  (1) 

f=-Vmcos(q  —  Om)  +  Vtcos(q-dt),  (2) 

em  =  am/Vm ,  (3) 

0t=at/Vt,  (4) 


where  r  and  r  are  the  relative  distance  and  the  relative  velocity 
between  missile  and  target,  q  and  q  are  the  LOS  angle  and  LOS  an¬ 
gular  rate  between  missile  and  target,  Vm  and  Vt  are  the  missile 
velocity  and  the  target  velocity,  6m  and  8t  are  the  flight-path  angle 
of  missile  and  target,  am  and  at  are  the  missile  lateral  accelera¬ 
tion  normal  to  missile  velocity  and  the  target  lateral  acceleration 
normal  to  target  velocity,  respectively. 

Differential  Eq.  (1)  with  respect  to  time,  yields 

rq  =  —2  rq  —  am  cos  (q  —  6m)  +  Vm  sin(q  —  0m) 

+  atcos(q -6»t)  -  VfSin(q -6>t).  (5) 

The  design  objects  here  are  to  develop  guidance  laws  such  that 
the  guidance  law  can  not  only  guarantee  the  missile  hits  the  tar¬ 
get  with  zero  miss  distance,  but  also  guarantee  the  missile  has  an 
appropriate  attitude  to  intercept  the  maneuvering  target.  That  is 
to  design  the  guidance  laws  am  such  that  the  LOS  angular  rate  q 
will  converge  to  zero,  and  the  LOS  angle  q  will  converge  to  a  de¬ 
sired  LOS  angle.  The  main  results  will  be  shown  in  the  following 
sections. 

In  this  paper,  the  signals  r,  r,  q,  q  and  0m  are  assumed  to  be 
measurable.  Let  qd  be  as  the  desired  LOS  angle.  Let  qe  be  the  LOS 
angle  error,  i.e.,  qe  =  q  —  qd.  Defining  xi  =  qe,  x2  =  qe,  Eq.  (5)  can 
be  written  as 

|*i  =x2, 

*2  =  7  (-2 rx2  -  am  cos(q  -  dm)  +  Vm  sin (q  -  6m )  (6) 

+  at  cos (q  -  0t )  -  Vt  sin(q  -  0t)). 

Assumption  1.  The  target  lateral  acceleration  at(t)  and  longitu¬ 
dinal  acceleration  Vt(t)  are  assumed  to  be  bounded  and  satisfy 
|at (t)|  <  d i,  |Vt(t)l  <  d2,  for  all  t  >  0,  where  di  and  d2  are  the 
known  bound  of  the  target  accelerations. 

Assumption  2.  During  the  time  horizon  of  the  guidance  process, 
it  has  f(t)  <  0,  0  <  r(t)  <  r(0),  t  >  0  [46],  The  missile  intercepting 
target  by  impact  (“hit-to-kill”)  happens  when  r  /  0  but  belongs  to 
the  interval  r°  e  [rmin,  rmax]  =  [0.1,  0.25]  m  [30]. 

Before  giving  the  guidance  laws  design,  some  results  on  the  fi¬ 
nite  time  stability  are  needed.  Firstly,  the  definition  of  finite-time 
stability  is  presented  as  follows. 

Definition  1.  (See  [3].)  Consider  a  nonlinear  system  in  the  form 
*=/(*),  /(0)  =  0,  X£R", 


Lemma  1.  (See  [3].)  Consider  the  nonlinear  system  described  by  Eq.  (7). 
Suppose  there  exists  a  continuous  function  V(x)  :U  -*■  R  such  that: 

(i)  V  (x)  is  positive  definite. 

(ii)  There  exist  real  numbers  c  >  0  and  a  £  (0, 1)  and  an  open  neigh¬ 
borhood  Uq  C  U  of  the  origin  such  that  V(x)  <  — cV“(x),  x  £ 
Uo\{0). 

Then,  the  origin  is  a  finite-time  stable  equilibrium  of  system  (7).  More¬ 
over,  if  T  is  the  settling  time,  then  T(x)  ^  V1~“(x)  for  all  x  in 

some  open  neighborhood  of  the  origin.  IfU  —  Uo  —  Rn,  the  origin  is  a 
globally  finite-time  stable  equilibrium  of  system  (7). 

3.  Guidance  law  based  on  integral  sliding  mode  control 

3.1.  Guidance  law  based  on  linear  integral  sliding  mode  control 

In  this  section,  a  linear  ISM  guidance  law  is  presented  for  the 
terminal  guidance  system  described  by  Eqs.  (l)-(4)  and  (6).  Let 
us  give  a  brief  description  of  the  design  method.  Firstly,  a  linear 
ISM  manifold  inspired  by  Zong  [49]  is  designed.  Secondly,  a  linear 
ISM  guidance  law  based  on  the  linear  ISM  manifold  is  developed, 
and  then  it  is  proved  that  the  obtained  guidance  law  can  guar¬ 
antee  that  the  LOS  angular  rate  asymptotically  converges  to  zero, 
and  guarantee  that  the  LOS  angle  asymptotically  converges  to  the 
desired  LOS  angle. 

Considering  the  process  of  the  terminal  guidance  of  missile,  the 
following  linear  ISM  manifold  is  introduced 


S]  =x2 -x2(0)+  /  (kiXi  +k2x2)dt,  0,  (8) 


where  k\  >  0,  k2  >  0  are  parameters  to  be  designed.  Note  that  it  is 
Si  (0)  =  0  at  t  =  0,  that  means  the  LOS  angle  error  qe  and  the  LOS 
angular  rate  error  qe  occur  on  the  ISM  manifold  (8)  starting  from 
the  initial  time  instance. 

Theorem  1.  Consider  guidance  system  and  .If 

1  and  2  are  satisfied,  and  the  switch  gain  satisfies  e  >  yjd\  +  d\,  the 
following  linear  ISM  guidance  law  based  on  the  linear  ISM  manifold  (8) 

_  — 2rx2  +  kirxi  +  k2rx2  +  Vm  sin(q  -0m)+e  sign(si) 


can  guarantee  that  the  LOS  angular  rate  q  asymptotically  converges  to 
zero,  and  the  LOS  angle  q  asymptotically  converges  to  the  desired  LOS 
angle  qd. 

Note,  sign(-)  in  linear  ISM  guidance  law  (9)  denotes  the  standard 
(7)  signum  function. 
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Proof  of  Theorem  1.  Differentiation  of  Eq.  (8)  with  respect  to  time 
yields 


si  =  x2  +  k2x2  +  kixi.  (10) 

Substitution  of  Eqs.  (6)  and  (9)  into  Eq.  (10)  yields 


Si  =  -  [-2 rx2  -  am  cos(q  -  dm)  +  Vm  sin(q  -  0m ) 

+  at  cos(q  -  0t)  -  Vt  sin(g  -  dt)]  +  k2x2  +  Mi 
=  1  [at  cos(q  -  9t)  -  Vt  sin(q  -Gt)-s  sign(s)].  (11) 

Consider  the  final  Lyapunov  function 

Vi  =  \s2v  (12) 

The  derivative  of  Eq.  (12)  along  the  trajectories  of  Eq.  (11).  yields 
l/i  =  si  si  =  y  [at  cos(q  -  0t)  -  Vt  sin(q  -  0t)  -  s  sign(si)] 


V  Qt  +  v? 

=  —  [Jat  +  V?  cos((3  -&t  +  4>)-£  Sign  (Si )] 

s-Jd\  +  d\ 


According  to  the  Routh  Criterion,  the  condition  fei  >  0,  k2  >  0  is 
sufficient  to  guarantee  that  the  LOS  angle  error  qe  asymptotically 
converges  to  zero  with  infinite  convergence  time.  Hence,  the  LOS 
angular  rate  error  qe  asymptotically  converges  to  zero.  Note  that 
qe  =  q  —  <jd  =  q  and  qe  =  q  —  qd.  which  demonstrate  that  the  LOS 
angular  rate  q  will  asymptotically  converge  to  zero  and  the  LOS 
angle  q  will  asymptotically  converge  to  the  desired  LOS  angle  qd- 
The  proof  is  finished.  □ 

3.2.  Guidance  law  based  on  nonlinear  integral  sliding  mode  control 

The  proposed  linear  ISM  guidance  law  in  Section  3.1  can  guar¬ 
antee  the  LOS  angular  rate  and  the  LOS  angle  asymptotical  con¬ 
vergence  with  infinite  settling  time.  To  further  improve  the  con¬ 
vergence  characteristics  of  LOS  angle  and  LOS  angular  rate,  one 
possible  method  is  to  introduce  a  nonlinear  control  strategy.  In 
this  section,  a  finite-time  control  strategy  is  introduced.  Firstly,  a 
nonlinear  ISM  manifold  based  on  a  homogeneity  controller  for  in¬ 
tegrator  system  is  introduced.  And  then,  a  nonlinear  ISM  guidance 
law  based  on  the  nonlinear  ISM  manifold  is  designed.  Finally,  it  is 
proved  that  the  obtained  guidance  law  can  ensure  LOS  angle  and 
LOS  angular  rate  the  finite-time  convergence  characteristics.  Be¬ 
fore  the  design  of  finite-time  convergence  guidance  laws,  we  need 
some  results  on  the  finite-time  stability  for  nonlinear  systems. 

Lemma  2.  (See  [3].)  For  a  e  (0, 1),  the  feedback  law 

u  =0(fi,?2)  =  — fci Ifi  1“  sign(^i)  -k2|?2|^r  sign (£2)  (19) 

renders  the  origin  finite-time  stable  for  the  double  integrator  system 

tl=f2,  ll  =  U,  (20) 


where  cos  d>  —  ,  a‘  and  sin</>  =  ,  v‘ 

According  to  Assumption  2,  during  the  time  horizon  of  the 
guidance  process,  it  has 


where  k\  >0,k2>  0. 

Remark  1.  (See  [14].)  The  origin  of  system  (20)  is  finite-time  stable. 
There  exists  an  explicit  finite-time  Lyapunov  function 


r(t)  <0,  0  <  r(t)  <  r(0),  t  >  0. 


(14) 


Combining  Eqs.  (13)  and  (14),  it  obtains 


Vi 


r(0)  1  11  r(0)  1  ' 


(15) 


According  to  Lemma  1,  Eq.  (15)  indicates  that  the  LOS  angle 
error  qe  converges  to  the  linear  ISM  manifold  (8)  in  finite-time. 
Assume  that  the  LOS  angle  error  qe  reaches  the  manifold  at  ti. 
According  to  Lemma  1,  it  can  obtain 


V2r(0) 


=  |si  (0)  |  =  0.  (16) 


Eq.  (16)  demonstrates  that,  even  if  there  exists  the  external  dis¬ 
turbance,  the  LOS  angle  error  qe  and  the  LOS  angular  rate  error 
qe  occur  on  the  ISM  manifold  (8)  starting  from  the  initial  time  in¬ 
stance. 

In  the  end,  we  prove  the  LOS  angle  error  qe  can  asymptotically 
converge  to  zero.  Because  of  the  system  states  qe  and  qe  occurring 
on  linear  ISM  manifold,  it  follows  from  Eq.  (8)  that 


Si  =x2 -x2(0)  +  y  (kiXi  +k2x2)dt  =  0,  (17) 

0 

which  implies  that  Si  =  0,  i.e., 


si  =  x2  +  k2x2  +  kiXi  =  qe  +  k2qe  +  k1qe  =  0.  (18) 


9  4-  9l43+a0/2 

Vo=  3+a  lfl|(3+g)/2 

+  ^^lf2|(3+“)/(1+“)+k0f1f2, 

5a  + 1 


(21) 


where  ko  =  (ki  /k2)1^“  >  0  with  k2  sufficiently  large,  and  then  it 
can  obtain  that  there  is  a  K  such  that 


V0  <  _k-^2+2“)/(3+“).  (22) 

From  Lemma  1,  this  inequality  leads  to  the  convergence  time  esti¬ 
mation  given  by 


T  ^  (3  +  a)V^a)/(3+a)($o) 

K(  1-a) 

where  £0  is  the  initial  state  of  the  system. 


(23) 


Similarly  to  the  design  method  of  linear  ISM  guidance  law, 
a  nonlinear  ISM  manifold  inspired  by  Zong  [49]  is  introduced, 

S2=x2-x2(0)  +  J  [ki  | x  1 1“  signfxi )  +  k2 |x2 1  sign(x2)]  dt, 

0 

0,  (24) 

where  ki  >  0,  k2  >  0,  0  <  a  <  1  are  parameters  to  be  designed. 
Note  that  it  is  also  s2(0)  =  0  at  t  =  0,  that  means  the  LOS  angle 
error  qe  and  the  LOS  angular  rate  error  qe  occur  on  the  nonlinear 
ISM  manifold  (24)  starting  from  the  initial  time  instance. 
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Theorem  2.  Consider  guidance  system  (l)-(4)  and  (6).  If  Assumptions 
1  and  2  are  satisfied,  and  the  switch  gain  satisfies  e  >  Jd\  +  d\,  the 
following  nonlinear  ISM  guidance  law  based  on  the  nonlinear  ISM  man¬ 
ifold  (24) 

am  =  (— 2rx2  +  k]  r|xi  |“  sign(xi)  +k2r\x2\<*+'  sign(x2) 

+  Vm  sin(q  -  9m)  +  esign(s2))(cos(q  -  6»m))_1  (25) 

can  guarantee  that  the  LOS  angular  rate  q  converges  to  zero  infinite  time, 
and  the  LOS  angle  q  converges  to  the  desired  LOS  angle  qd  infinite  time. 


Proof.  Differentiation  of  Eq.  (24)  with  respect  to  time  yields 


s2  =  x2+k2 |x2| sign(x2)+ki|xi|“sign(xi).  (26) 

Substitution  of  Eqs.  (6)  and  (21)  into  Eq.  (22)  yields 
s2  =  ^  [— 2rx2  -  am  cos(q  -  9m )  +  Vm  sin(q  -  9m) 

+  at  cos(q  —  9t)  —  Vt  sin(q  —  0t)] 

+  k2  |x2 1  sign(x2)  +  ki  |xi  sign(xi ) 

=  ^[at  cos (q  -  9t)  -  Vt  sin (q  -  9t )  -  esign(s2)].  (27) 

Similarly  to  the  proof  of  Theorem  1,  the  final  Lyapunov  function  is 
chosen  as  V2  =  \s^.  The  derivative  of  Lyapunov  function  along  the 
trajectory  of  Eq.  (27),  yields 


V2  =  s2s2  ^ 


V2(e 


v\'2. 


(28) 


According  to  Lemma  1,  Eq.  (28)  demonstrates  that  the  LOS  an¬ 
gle  error  qe  converges  to  the  nonlinear  ISM  manifold  (24)  in  finite 
time.  Assume  that  the  LOS  angle  error  reaches  the  manifold  at  t2. 
According  to  Lemma  1,  it  can  obtain 


t2  < - r,(0)  |s2(0)|  =  0.  (29) 

e-Jd2,+d22 

Eq.  (29)  demonstrates  that,  even  if  there  exists  the  external  distur¬ 
bance,  the  LOS  angle  error  qe  and  the  LOS  angular  rate  error  qe 
occur  on  the  nonlinear  ISM  manifold  (24)  starting  from  the  initial 
time  instance. 

In  the  end,  we  prove  that  the  LOS  angle  error  qe  will  converge 
to  zero  in  finite  time.  Because  of  the  system  states  qe  and  qe  oc¬ 
curring  on  nonlinear  ISM  manifold,  it  follows  from  Eq.  (24)  that 


According  to  Remark  1,  the  convergence  time  estimation  is 
given  by 


'  K(l-o) 


where  K  >  0,  Xo  is  the  initial  state  of  the  system. 
The  proof  is  finished.  □ 


(32) 


Remark  2.  On  the  one  hand,  due  to  existence  of  function  esign(s), 
both  ISM  guidance  laws  (9)  and  (25)  are  discontinuous,  thus  the 
system  may  generate  chattering  problem.  On  the  other  hand,  to 
guarantee  that  the  guidance  system  has  a  good  disturbance  rejec¬ 
tion  performance,  the  switch  gain  e  needs  to  be  chosen  larger  than 
the  bound  of  the  target  acceleration.  However,  if  there  is  no  precise 
information  on  the  bound  of  the  target  acceleration,  the  switching 
gain  will  be  taken  to  be  large  enough.  This  inaccuracy  worsens 
the  chattering  caused  by  sliding  mode  control.  Note  that,  the  func¬ 
tions  fcir|xi  |  sign(xi)  and  k2r|x2|“  sign(x2)  in  guidance  law  (25) 
are  continuous,  which  do  not  generate  chattering.  To  handle  chat¬ 
tering  phenomenon,  one  of  the  solution  is  to  employ  a  saturation 
function  es/(|s|  +  S),  where  &  is  a  very  small  positive  number,  to 
substitute  the  signum  function  esign(s)  in  guidance  laws  (9)  and 
(25).  This  method  can  remove  the  chattering  phenomenon.  How¬ 
ever,  to  do  this,  the  disturbance  rejection  performance  is  sacrificed 
to  some  extent.  Another  efficient  method  for  suppressing  the  chat¬ 
tering  phenomenon  is  to  employ  a  NDOB  to  estimate  the  target 
acceleration.  This  estimate  of  target  acceleration  will  be  employed 
as  a  compensation  item  in  ISM  guidance  laws.  After  disturbance 
compensation  based  on  NDOB,  the  switching  gain  only  needs  to 
be  larger  than  the  bound  of  the  disturbance  compensation  error 
which  is  usually  much  smaller  than  that  of  the  target  accelera¬ 
tion.  Hence,  the  chattering  phenomenon  will  be  reduced  while  the 
disturbance  rejection  performance  of  the  guidance  system  can  be 
maintained. 


4.  Composite  guidance  law  based  on  nonlinear  disturbance 
observer 

To  handle  the  chattering  phenomenon  and  maintain  the  perfor¬ 
mance  of  the  guidance  system,  we  will  integrate  the  ISM  guidance 
laws  and  NDOB  to  give  composite  guidance  control  schemes.  First 
of  all,  by  designing  a  non-smooth  NDOB,  the  target  acceleration 
can  be  estimated  in  finite  time.  This  estimated  target  accelera¬ 
tion  will  be  employed  as  a  feedforward  compensation  item  in  the 
ISM  guidance  laws  (9)  and  (25).  Then,  composite  guidance  control 
schemes  can  be  obtained. 


s2=x2-x2(0) 

+  /fc*r  sign(xi )  +  k2 |x2 1 sft  sign(x2)]  dt  =  0, 
o 

t  >  t2  =  0,  (30) 

which  implies  that  s2  =  0,  i.e., 
s2  =  x2  +  k2 |x2|^pt  sign(x2)  +  ki  f*i  |“  sign(xi) 

=  qe  +  k2  Iqel^Ff  sign(qB)  +  ki|qe|“sign(qe)  =  0,  (31) 

where  k\  >  0,  k2  >  0,  0  <  a  <  1. 

According  to  Lemma  2,  system  (31)  is  finite-time  stable.  There¬ 
fore,  the  LOS  angle  error  qe  and  the  LOS  angular  rate  error  qe 
converge  to  zero  in  finite  time.  Note  that  qe  =  q  —  qd  =  q  and 
qe=q  —  qd ,  which  mean  the  LOS  angular  rate  q  will  converge  to 
zero  in  finite  time  and  the  LOS  angle  q  will  converge  to  the  desired 
LOS  angle  qd  in  finite  time. 


4.1.  Nonlinear  disturbance  observer 

Consider  single-input  single-output  (SISO)  dynamics 
cr=g(t)  +  U,  (33) 

where  cr  e  R,  u  e  R  is  a  control  input,  g(t)  is  an  uncertain  suf¬ 
ficiently  smooth  function  and  satisfies  g(t)  ^  L,  L  >  0  is  Lip- 
shitz  constant.  The  control  function  u(t)  is  Lebesgue  measurable. 
Eq.  (33)  is  understood  in  the  Filippov  sense  [18],  which  means  in 
particular  that  cr  (t)  is  an  absolutely  continuous  function  defined 
for  any  t>  0.  By  a  simple  modification  of  Theorem  6  in  [18],  we 
have  the  following  result. 

Lemma  3.  (See  [18].)  Consider  the  following  NDOB: 

Iz0  =  v0  +  u,  v0  =  — A.0|z0  -  a  |2/3  sign(z0  -  cr)  +  z\ , 
zi  =  vi,  vi  =  — A.i  |zi  -  v0|1/2sign(zi  -  v0)  +  z2,  (34) 

z2  =  -X2  sign(z2  — vi), 
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where,  A,  >  0,  i  —  0, 1 , 2  are  sufficiently  large,  then  it  can  obtain  that  Z\ 
converges  to  g(t)  in  finite  time. 

Remark  3.  It  is  proved  in  [18]  that  in  order  to  obtain  the  finite¬ 
time  convergence  of  the  observation  errors,  the  parameters  A,-, 
i  =  0,1,2  should  be  chosen  large  enough.  However,  according 
to  [29],  the  simulation  shows  that  Ao  =  2 11/3,  Ai  =  1.5L1/I2  and 
A3  =  1.1L  are  enough  to  ensure  the  stability  of  observation  error. 


am  =  (-2  rx2  +  kirxi  +  k2rx2  +  Vm  sin  (q  -6m)  +  r}  sign(si) 

+  Wq(t))(cOS(q  —  dm))1  (38) 

can  guarantee  that  the  LOS  angular  rate  q  asymptotically  converges  to 
zero,  and  the  LOS  angle  q  asymptotically  converges  to  the  desired  LOS 
angle  qd. 

Proof.  Substituting  Eqs.  (6)  and  (38)  into  Eq.  (10),  yields 


According  to  guidance  system  (6),  define  wq  =  atcos(q  —  Of)  — 
Vt  sin(q  -  9t)  which  is  the  target  acceleration  component  normal 
to  the  LOS.  The  design  object  is  to  develop  a  NDOB  such  that  the 
estimated  target  acceleration  wq  can  converge  to  actual  target  ac¬ 
celeration  Wq  in  finite  time. 

And,  define  d(t)  =  ywq  as  lumped  disturbance  and  d(t)  as  the 
estimated  disturbance.  Let  G(t)  =  g(t)  +  d(t),  g(t)  =  l(-2 rx2  + 
Vm  sin(q  —  6m))  and  B  =  —  J  cos(q  —  6m).  Guidance  system  (6)  can 
be  written  as 


=x2, 

{  x2  =  G(t)  +  Bam. 


(35) 


The  derivatives  of  functions  g(t)  and  d(t)  are  g(t)  —  \[2 r2x2  — 
2 rrx2  -  2 rrx2  +  (rVm  -  rVm)  sin(q  -  6m)  +  rVm(q  -  0m)  cos(q  -  6>m)], 
d(t)  =  p  [(rdt  -  rat)  cos(q  -  9t)  +  (rVt  -  rVt)  sin(q  -  0t)  -  r(q  -  9t)  x 
(Vt  cos(q  -  0t)  +  at  sin(q  -  (9t))],  respectively. 


si  =  -  [ w,  (t)  -  w,  (t)  -  r\  sign(si )] .  (39) 

Similarly  to  the  proof  of  Theorem  1,  the  final  Lyapunov  function 
is  chosen  as  V3  =  ^s2.  The  derivative  of  Lyapunov  function  along 
the  trajectory  of  Eq.  (39),  yields 


V3  =  SiSl  =  —  [Wq(t)  -  Wq(t)  -  T)  Sign(Si)]  - - -  |Si  | 


V2Ql  -  Q)  y\/2 
r(0)  3 


(40) 


According  to  Lemma  1,  Eq.  (40)  demonstrates  that  the  LOS  angle 
error  qe  converges  to  the  linear  ISM  manifold  Eq.  (8)  in  finite  time. 

The  rest  procedure  is  identical  with  that  of  Theorem  1 ,  which  is 
omitted  here.  Similarly  to  the  proof  of  Theorem  1,  we  obtain  that 
the  LOS  angular  rate  q  will  asymptotically  converge  to  zero  and  the 
LOS  angle  q  will  asymptotically  converge  to  the  desired  LOS  angle 
qd.  The  proof  is  finished.  □' 


Remark  4.  If  the  target  acceleration  at  and  Vt  are  differentiable, 
then  the  function  g(t)  and  d(t)  are  differentiable  except  for  r  =  0. 
According  to  Assumption  2,  the  functions  g(t)  and  d(t)  are  dif¬ 
ferentiable  and  the  g(t)  and  d(t)  are  continuous  everywhere  until 
hit-to-kill  happens.  Since  the  function  G(t)  =  g(t)  +  d(t)  is  continu¬ 
ously  differentiable  for  any  r  ^  r°,  its  derivative  G(t)  has  a  Lipshitz 
constant  L. 

For  guidance  system  (6),  according  to  Lemma  3,  the  lump  dis¬ 
turbance  d(t)  can  be  estimated  by  the  following  NDOB  (36) 

z0  =  v0  +  Bam  +  g(t), 

Vo  =  -Aolzo  -  x2|2/3  sign(z0  -  x2)  +  zi , 
zi=vi,  Vl  =  — Ai  |zi  -  vo|1/2sign(z1  -  vQ)  +  z2,  (36) 

z2  =  -A2  sign(z2  —  Vi), 
d(t)  =  zi . 

According  to  Lemma  3  and  Remark  3,  it  can  obtain  that  the  es¬ 
timated  lumped  disturbance  d(t)  converges  to  d(t)  in  finite  time, 
i.e.,  the  estimated  target  acceleration  wq(t)  =  rd(t )  converges  to 
real  target  acceleration  wq(t)  in  finite  time. 

4.2.  Composite  guidance  law 


Theorem  4.  Consider  guidance  system  and  .  If 

2  and  3  are  satisfied,  and  the  switch  gain  satisfies  q  >  q,  the  following 
nonlinear  composite  guidance  law  combining  the  nonlinear  ISM  guid¬ 
ance  law  (25)  with  NDOB  (36) 

am  =  (-2 rx2  +  kir|xi  |“  sign(xi)  +  k2r|x2| sign(x2) 

+  Vm  sin(q  -  0m)  +  q  sign(s2)  +  w,(t)) (cos(q  -  6>m))_1 

(41) 

can  guarantee  that  the  LOS  angular  rate  q  converges  to  the  origin  in  finite 
time,  and  the  LOS  angle  q  converges  to  the  desired  LOS  angle  qd  infinite 


Proof.  Substituting  Eqs.  (6)  and  (41)  into  Eq.  (26),  yields 

52  =  ^[wq(t)  -  Wq(t)  -  q  sign(s2)] .  (42) 

Similarly  with  the  proof  of  Theorem  2,  the  final  Lyapunov  func¬ 
tion  is  chosen  as  V4  =  ^s2,  also.  The  derivative  of  Lyapunov  func¬ 
tion  along  the  trajectory  of  Eq.  (42),  yields 


V2(>?  ~  Q)  „i/2 


(43) 


Before  giving  the  composite  guidance  law,  the  following  as¬ 
sumption  is  presented  on  the  basis  of  Lemma  3. 

Assumption  3.  The  target  acceleration  estimation  error  is  bounded 
and  there  exists  a  constant  Q,  such  that 

|Wq(t)-Wq(t)|<£.  (37) 


According  to  Lemma  1,  Eq.  (43)  demonstrates  that  the  LOS  angle 
error  qe  converges  to  the  nonlinear  ISM  manifold  Eq.  (24)  in  finite 

The  rest  procedure  is  identical  with  that  of  Theorem  2,  which 
is  omitted  here.  Similarly  to  the  proof  of  Theorem  2,  we  obtain 
that  the  LOS  angular  rate  q  will  converge  to  zero  in  finite  time  and 
the  LOS  angle  q  will  converge  to  the  desired  LOS  angle  qd  in  finite 
time.  The  proof  is  finished.  |I: 


Theorem  3.  Consider  guidance  system  (1 )— (4)  and  (6).  If  Assumptions 
2  and  3  are  satisfied,  and  the  switch  gain  satisfies  q  >  Q,  the  following 
composite  linear  guidance  law  combining  the  linear  ISM  guidance  law 
(9)  with  NDOB  (36) 


Remark  5.  The  differences  between  Theorems  1,  2  and  Theorems 
3,  4  are  that  feedforward  compensation  terms  wq  are  introduced 
into  controller  in  Theorems  3  and  4.  When  the  parameters  of  the 
NDOB  are  appropriately  chosen,  wq  will  be  well  estimated.  Usually 
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the  bound  of  disturbance  compensation  errors  can  be  less  than  the 
disturbances.  In  this  case,  the  composite  guidance  laws  (38)  and 
(41)  may  select  a  smaller  value  for  the  switching  gain  r)  without 
sacrificing  disturbance  rejection  performance,  which  helps  to  re¬ 
duce  large  chattering  caused  by  high  gain. 


The  response  of  LOS  angle,  LOS  angular  rate,  estimated  distur¬ 
bance,  missile  acceleration  command,  sliding  mode  manifold  and 
interception  geometry  are  shown  in  Figs.  2(a),  (b),  (c),  (d),  (e) 
and  (f),  respectively.  The  miss  distances  and  interception  time  are 
shown  in  Table  1. 


Remark  6.  The  composite  guidance  law  and  NDOB  are  designed 
separately.  The  control  parameters  and  the  NDOB  parameters  can 
be  also  chosen  separately.  The  NDOB  parameters  k,-,  i  —  0, 1, 2  can 
be  obtained  according  to  Remark  3.  For  the  Lipshitz  constant  L, 
when  L  is  chosen  from  small  to  large  value,  the  disturbance  es¬ 
timation  dynamics  response  becomes  faster.  Flowever,  when  its 
value  is  large  enough,  the  overshoot  appears  accordingly.  The  pa¬ 
rameters  L  should  be  selected  suitably  considering  trade  off  be¬ 
tween  response  speed  and  overshoot.  The  guidance  law  parameters 
fci  >  0,  k2  >  0,  0  <  a  <  1,  s  >  0  and  rj  >  0  directly  influence  the 
closed  loop  performance.  The  control  parameters  k \ ,  k-2  can  be  ob¬ 
tained  by  pole  assignment  method.  Form  the  stability  analysis  of 
Theorems  1  to  4,  switch  gain  e  should  be  larger  than  the  bound  of 
lumped  disturbances  for  ISM  guidance  laws  (9),  (25),  and  switch 
gain  rj  should  be  large  than  the  bound  of  disturbance  compensa¬ 
tion  error  for  composite  guidance  laws  (38),  (41).  Usually,  when 
the  parameter  a  is  chosen  smaller,  the  convergence  rate  of  states 
will  be  faster. 


5.  Simulation  results 

In  this  section,  the  numerical  examples  are  performed  to  il¬ 
lustrate  the  performance  of  all  proposed  guidance  laws  (9),  (25), 
(38),  (41).  The  initial  states  of  guidance  system  are  chosen  as  fol¬ 
lows.  The  initial  relative  distance  between  missile  and  target  is 
ro  =  5000  m.  The  initial  LOS  angle  is  qo  =  30  deg.  The  desired 
LOS  angle  is  qd  =  20  deg.  The  initial  missile  velocity  is  Vm0  — 
600  m/s  with  initial  flight-path  angle  9m o  =  60  deg.  In  the  sim¬ 
ulation,  the  missile  longitudinal  acceleration  is  Vm  —  0  m/s2,  i.e., 
the  missile  velocity  Vm  is  a  constant.  The  initial  target  velocity  is 
Vto  —  300  m/s  with  flight-path  angle  9to  =  0  deg.  The  accelera¬ 
tion  of  gravity  is  g  —  9.8  m/s2.  The  initial  states  of  NDOB  (36)  are 
chosen  as  Zo  =  [0, 0.5, 0]T. 

In  order  to  demonstrate  the  effectiveness  of  the  proposed  guid¬ 
ance  laws,  the  PN  guidance  law  is  chosen  for  comparison.  The  PN 
guidance  law  is  taken  as 

nc  =  -Nrq,  (44) 

where  N  is  a  unitless  coefficient  of  PN  guidance  law.  To  have  a  fair 
comparison  for  the  simulation  results,  the  acceleration  command 
is  limited  not  to  exceed  50  g.  The  parameter  N  in  PN  guidance 
law  is  chosen  as  N  =  5.  According  to  Remark  6,  the  parameters  ki 
and  fc2  of  all  proposed  guidance  laws  are  obtained  by  pole  place¬ 
ment  method.  The  parameters  of  the  linear  ISM  guidance  law  (9) 
and  composite  linear  guidance  law  (38)  are  selected  as  k\  =  1  and 
It2  =  2.  And,  the  parameters  of  the  nonlinear  ISM  guidance  law 
(25)  and  composite  nonlinear  guidance  law  (41)  are  selected  as 
ki  =  0.5,  k2  =  1  and  a  —  0.3.  The  switch  gain  e  of  ISM  guidance 
laws  (9)  and  (25)  is  selected  as  e  =  70.  And,  the  switch  gain  q 
of  composite  guidance  laws  (38)  and  (41)  is  selected  as  q  =  0.1. 
The  Lipshitz  constant  L  of  NDOB  is  chosen  as  L  =  50.  According 
to  Remark  3,  the  gains  Xo,  M  and  a 2  of  NDOB  can  be  chosen  as 
k0  =  7.368,  ki  =  10.607,  a2  =  55.  The  following  cases  are  simu¬ 
lated. 


Case  1.  It  is  supposed  that  the  target  acceleration  is  chosen  a 


at(t)  = 


=  (5g. 


(45) 


Note.  The  abbreviations  in  Fig.  2  and  Table  1,  i.e.,  PNGL,  LISMGL, 
NISMGL,  CLGL  and  CNGL,  denote  the  PN  guidance  law  (44),  linear 
ISM  guidance  law  (9),  nonlinear  ISM  guidance  law  (25),  composite 
linear  guidance  law  (38)  and  composite  nonlinear  guidance  law 
(41),  respectively.  In  the  following  cases,  the  abbreviations  have 
the  same  meanings. 

From  Table  1,  it  is  clear  that  four  proposed  guidance  laws  (9), 
(25),  (38),  (41)  can  guarantee  the  miss  distance  of  guidance  system 
less  than  0.01  m  in  spite  of  the  target  acceleration  with  signifi¬ 
cant  changes.  According  to  Ref.  [30],  it  indicates  that  the  missile 
can  intercept  the  large  maneuvering  target  by  hit-to-kill  guidance 
strategy.  Flowever,  the  PN  guidance  law  makes  the  missile  has  a 
large  miss  distance  353.64  m,  which  means  the  missile  will  miss 
the  target.  Therefore,  the  four  proposed  guidance  laws  are  available 
to  intercept  the  large  maneuvering  target,  while  the  PN  guidance 
law  cannot. 

As  shown  in  Fig.  2(a),  it  can  be  seen  that  the  four  proposed 
guidance  laws  can  guarantee  the  LOS  angle  converges  to  the  de¬ 
sired  angle  20  deg.  Fig.  2(b)  shows  that  all  proposed  guidance 
laws  can  guarantee  the  LOS  angular  rate  converges  to  zero.  But, 
the  PN  guidance  law  cannot  guarantee  the  convergence  of  LOS  an¬ 
gle  and  LOS  angular  rate.  From  Figs.  2(a')  and  2(b'),  it  is  clear 
that  the  nonlinear  guidance  laws  (25),  (41)  guarantee  LOS  an¬ 
gle  and  LOS  angular  rate  more  rapidly  convergence  rate  than 
that  of  linear  guidance  laws  (9),  (38).  And,  the  nonlinear  guid¬ 
ance  laws  (25),  (41)  can  guarantee  the  LOS  angle  and  LOS  an¬ 
gular  rate  finite-time  convergence.  Fig.  2(c)  shows  that  the  de¬ 
signed  NDOB  can  observe  the  target  acceleration.  It  is  clear  from 
Fig.  2(c')  that  the  estimated  target  acceleration  wq  converges  to 
the  real  target  acceleration  wq  in  finite  time.  The  control  input 
signals  of  all  proposed  guidance  laws  are  shown  in  Figs.  2(d,  d'). 
It  is  clear  that,  for  ISM  guidance  laws  (9),  (25),  the  guidance  com¬ 
mands  cause  the  chattering  phenomenon.  That  is  because,  to  sup¬ 
press  the  target  acceleration,  the  switch  gain  e  is  chosen  larger 
than  the  bound  of  target  acceleration.  But,  for  composite  guid¬ 
ance  laws  (38),  (41),  the  undesired  chattering  is  reduced  effec¬ 
tively.  That  is  because  that,  under  the  influence  of  NDOB,  the 
estimated  target  acceleration  wq  can  be  a  good  estimate  of  ac¬ 
tual  target  and  the  upper  bound  of  wq  —  wq  can  be  less  than 
Wq.  In  this  case,  the  switching  gain  q  required  usually  is  much 
smaller  than  e.  Thus,  the  discontinuous  terms  of  the  composite 
guidance  laws  (9),  (25)  can  be  much  smaller  and  the  chattering 
will  be  reduced.  Figs.  2(e,  e')  depict  the  performance  of  sliding 
mode  manifold.  It  can  be  seen  that  the  sliding  mode  is  stable  in 
spite  of  the  target  acceleration.  The  interception  geometries  are 
shown  in  Figs.  2(f,  f ),  which  shows  that  four  proposed  guidance 
laws  can  accomplish  the  interception,  while  the  PN  guidance  law 
cannot. 

In  order  to  demonstrate  that  the  control  strategy  can  also  work 
well  when  the  autopilot  dynamics  is  considered,  the  following  case 
is  simulated. 

Case  2.  On  the  basis  of  Case  1,  the  autopilot  dynamics  is  consid¬ 
ered  in  this  case.  During  the  design  process  of  guidance  laws,  the 
autopilot  dynamics  is  regarded  as  the  ideal  dynamics.  Flowever,  the 
autopilot  dynamics  exists  in  the  actual  missile  system.  Therefore, 
it  is  necessary  to  observe  the  interference  caused  by  autopilot  dy¬ 
namics.  In  this  case,  a  second-order  transfer  function  of  missile 
autopilot  dynamics  [17]  is  selected  as 
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Fig.  2.  Responses  under  five  guidance  laws  of  Case  1.  (a)  LOS  angle,  (b)  LOS  angular  rate,  (c)  Real  and  estimated  disturbances,  (d)  Missile  acceleration  command,  (e)  Sliding 
mode  manifold,  (f)  Interception  geometry. 


dm  S2  +  2  £C0nS  +Q)n 

where,  amo  is  the  actual  output  acceleration  of  missile  autopilot, 
f  is  the  damping  coefficient,  and  co„  is  the  natural  frequency,  re¬ 
spectively.  In  the  simulation  precess,  the  parameters  are  chosen  as 
co„  =  10,  f  =  0.5,  respectively.  The  response  of  LOS  angle,  LOS  an¬ 
gular  rate,  estimated  disturbance,  missile  acceleration  command, 
sliding  mode  dynamic  and  interception  geometry  are  shown  in 
Figs.  3(a),  (b),  (c),  (d),  (e)  and  (f),  respectively.  The  miss  distances 
and  interception  time  are  shown  in  Table  2. 

Similarly  to  Case  1,  it  can  be  observed  from  Table  2  that  four 
proposed  guidance  laws  (9),  (25),  (38),  (41)  can  guarantee  the  miss 
distance  of  guidance  system  less  than  0.01  m  while  the  autopilot 
dynamics  is  considered,  which  means  that  the  missile  can  accom¬ 
plish  hit-to-kill  guidance  strategy.  The  PN  guidance  law  makes  the 


Miss  distances  and  interception  time  of  Case  1. 

Guidance  law  Miss  distance  (m)  Interception  time  (s) 

PNGL  353.64  18.765 

LISMGL  6.79  xl0“12  17.077 

NISMGL  1.65  x  10~12  17.293 

CLGL  3.82  x  10~8  17.071 

CNGL  5.07  xl0“9  17.191 


missile  has  a  large  miss  distance  379.37  m,  which  means  the  mis¬ 
sile  will  miss  the  target.  Therefore,  the  four  proposed  guidance 
laws  are  available  to  intercept  large  maneuvering  target  in  the  ex¬ 
istence  of  autopilot  dynamics. 

Fig.  3(a)  depicts  that  the  four  proposed  guidance  laws  can  guar¬ 
antee  that  the  LOS  angle  converges  to  the  desired  angle  20  deg.  As 
shown  in  Fig.  3(b),  it  is  clear  that  proposed  guidance  laws  (25), 
(41)  can  guarantee  that  the  LOS  angular  rates  converge  to  zero. 
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Table  2 

Miss  distances  and  interception  time  of  Case  2. 


Guidance  law 


PNGL 

LISMGL 

NISMGL 

CLGL 

CNGL 


Miss  distance  (m) 


Interception  time 


379.37  18.787 

2.36  x  1(T2  17.287 

1.04  xl0“2  17.225 

2.12  x  1CT2  17.057 

4.18  x  10~2  17.116 


(s) 


And,  it  is  obvious  that  the  nonlinear  guidance  laws  guarantee  LOS 
angle  and  LOS  angular  rate  have  the  finite-time  convergence  rate, 
which  is  more  rapidly  than  that  of  linear  guidance  laws.  Figs.  3(c, 
d)  show  that  the  estimated  target  acceleration  converges  to  the 
real  target  acceleration  in  finite  time.  The  control  input  signals 
of  four  proposed  guidance  laws  are  shown  in  Figs.  3(d,  d').  By 
comparing  Figs.  2(d,  d')  and  Figs.  3(d,  d'),  one  can  observe  that, 
when  the  autopilot  dynamics  is  not  considered,  the  ISM  guidance 
laws  (9),  (25)  cause  heavily  chattering  phenomenon.  But,  when 


the  autopilot  dynamics  is  considered,  the  chattering  phenomenon 
will  be  moderated.  That  is  because  that  the  second-order  transfer 
function  acts  as  a  low-pass  filter  on  the  command.  It  will  re¬ 
duce  the  chattering  phenomenon  caused  by  discontinuous  term. 
It  is  clear  from  Fig.  3(e),  that  the  sliding  mode  manifolds  are 
affected  by  the  autopilot  dynamics.  The  interception  geometries 
are  shown  in  Fig.  3(f),  which  shows  that  four  proposed  guidance 
laws  accomplish  the  interception  while  the  PN  guidance  law  can¬ 
not. 

In  conclusion,  for  missiles  intercepting  large  maneuvering  tar¬ 
gets  with  impact  angle  constraint,  four  proposed  guidance  laws 
(9),  (25),  (38),  (41)  work  effectively  and  are  robust  to  target 
acceleration  and  autopilot  dynamics.  The  proposed  linear  guid¬ 
ance  laws  (9),  (38)  guarantee  the  LOS  angle  q  and  LOS  angular 
rate  q  asymptotical  convergence  property.  The  proposed  nonlin¬ 
ear  guidance  laws  (25),  (41)  employ  the  finite-time  control  strat¬ 
egy,  which  guarantee  the  LOS  angle  q  and  LOS  angular  rate  q 
finite-time  convergence  property.  Four  proposed  guidance  laws 


Z.  Zhang  et  al.  /  Aerospace  Science  and  Technology  3 1  (2013 )  30-41 


Fig.  3.  Responses  under  five  guidance  laws  of  Case  2.  (a)  LOS  angle,  (b)  LOS  angular  rate,  (c)  Real  and  estimated  disturbances,  (d)  Missile  acceleration  command,  (e)  Sliding 
mode  manifold,  (f)  Interception  geometry. 


can  ensure  that  the  missiles  hit  the  target  with  zero  miss  dis¬ 
turbance.  The  NDOB  can  work  effectively  to  estimate  the  target 
acceleration.  The  estimated  target  acceleration  converges  to  the 
real  target  acceleration  in  finite  time.  After  disturbance  compen¬ 
sation  based  on  NDOB,  composite  guidance  laws  (38),  (41)  have 
the  ability  to  reduce  the  chattering  and  remain  the  performance 
of  that  which  ISM  guidance  laws  (9),  (25)  have.  Thus,  the  pro¬ 
posed  guidance  laws  (9),  (25),  (38),  (41)  achieve  the  design  ob¬ 
jective. 

6.  Conclusion 

In  this  paper,  based  on  ISMC  and  NDOB  theory,  guidance  laws 
with  impact  angle  constraint  have  been  developed.  All  guidance 
laws  are  proposed  to  achieve  the  hit-to-kill  guidance  strategy 
in  the  presence  of  target  maneuvers.  Proposed  linear  guidance 
laws  can  make  the  LOS  angle  and  LOS  angular  rate  asymptot¬ 


ical  convergence.  Proposed  nonlinear  guidance  laws  can  guar¬ 
antee  the  LOS  angle  and  LOS  angular  rate  finite-time  conver¬ 
gence.  With  the  help  of  NDOB  by  estimating  the  target  ac¬ 
celeration,  the  composite  guidance  laws  are  proposed  to  re¬ 
duce  the  chattering  phenomenon  which  exists  in  ISM  guidance 
laws  and  to  remain  the  performance  of  that  which  ISM  guid¬ 
ance  laws  have.  Finally,  simulation  comparison  results  are  pro¬ 
vided  to  demonstrate  the  effectiveness  of  the  presented  meth¬ 
ods. 
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